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Abstract 

We re-examine the neutrino decay solution to the solar neutrino problem in the 
light of the SuperKamiokande (SK) data. For the decay solution the SK spectrum 
data by its own can provide a fit comparable to the fit obtained from the MSW 
solution. However when one combines the results from the total rates of the ^"^Cl 
and ''^Ga experiments the fit becomes much poorer. 
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In this paper we examine the status of the neutrino decay solution to the solar 
neutrino problem in the light of the SK data []. Neutrino decay as a solution to the 
solar neutrino problem has been considered earlier (pre-SK) 0, |^. However since 
the Ga data implies that the low energy pp neutrinos are less suppressed compared 
to the high energy suppression seen in CI or the Kamiokande experiments, this 
solution was ruled out at 99% C.L. since the decay term exp{—aL/E), (where 
a = m2/To, m2 being the mass of the unstable state and Tq its rest frame lifetime), 
suppresses the low energy neutrino flux more than the high energy flux. However 
the SK spectrum data shows more events for the high energy bins. Though the 
statistics in these high energy bins need improvement, we explore the status of the 
neutrino decay scenario in the context of this behaviour of the SK spectrum data. 

Radiative decays of neutrinos are severely constrained and thus we are in- 
terested in the non-radiative decay modes. Two classes of models have been con- 
sidered in the literature in this connection. 

1. In neutrinos are considered to be Dirac particles. They consider the decay 
mode z/2 uir + 0, where uir is a right handed singlet and (p is an iso-singlet 
scalar. Thus all the final state particles are sterile. 

2. In 1^ neutrinos are assumed to be Majorana particles and the decay mode 
is 1^2 ^ i^i + J, where z/i interacts as a z/g with a probability |?7eiP and J is 
a Majoron. 

We work with just two flavors for simplicity and assume that the state 1^2 is 
unstable, which decays with a rest frame lifetime Tq. The other neutrino mass 
states have lifetimes much greater than the sun-earth transit time and hence can 
be taken as stable. In presence of decay, 

P.... = (l-|f/e2r)' + |f/e2|'exp(-47rL/Arf) 

+ 2|f/e2|'(l - \Ue2\^) exp(-27rL/Arf) cos(27rL/Ao,e) (1) 



Pu.ux = |?7e2|'(l- |f/e2n{l + exp(-47rL/A,) 

- 2 exp(-27rL/Ad) cos(27rL/Ao,,)} (2) 

where x can be either /i or r in the two flavour framework. A^ is the decay length 
defined as, 

Xd = 2.5 X IQ-^km 3 

MeV a ^ ' 



^The possibility of neutrino decay as a solution to the atmospheric neutrino problem has been 
considered recently in |j 
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Xosc is the oscillation wavelength defined as 

A„,, = 2.5 X lO-^fcm-^^ (4) 
MeV Am^ 

L = -R(t)(l — r/R(t)), r being the distance of the point of neutrino production 
from the center of the sun and R{t) is the sun-earth distance given by, 



R{t) = Ro 



, t , 
1 - ecos(27r-) 



(5) 



Here, i?o = 1-49 x 10 cm is the mean sun-earth distance and e = 0.0167 is the 
ellipticity of the earth's orbit, t is the time of the year and T is 1 year. The Am^ 
dependent oscillations are important around 10"^'' — 10~^^ eV^. We assume Am^ 
to be much higher so that the cosine term averages out to zero. We further assume 
t/ei = cos Q and [/e2 = sin Q so that the survival probability is 

P,,,, = cos^ B + sin^ B exp(-47rL/Ad) (6) 

In Fig. 1 we plot Py^y^ as a function of a for sin^ Q = 0.51 for two illustrative 
values of z/g energies 7 and 13 MeV. For a < lO^^^eV^^, the exp(— 47rL/Ad) ~ 1 
signifying very little decay. Beyond a = lO^^^eV"^ the exp{—4nL/Xd) ~ which 
is the case where we have very fast decay. In both these regions the probabilities 
are energy independent. In the region where a is around lO'^^eV"^ the probability 
does depend on the energy. 

The details of the solar neutrino code employed for performing the x^-analysis 
is described in J^, ^. We use the BP98 solar model ||^ as the reference SSM. We 
perform x^-analyses 

• using the total rates from ^^C/, '^^Ga and SK experiments. We incorporate 
the theory errors and their correlations. 



• using the 825 day SK spectrum data [jT2[ including the uncorrelated as well 
as the energy-bin correlated errors. 

• of the combined rates and SK spectrum data 

For the last two cases we vary Xb - the normalization of the fiux with respect 
to the SSM value as free parameter and determine its best-fit value. 

From the fact that the neutrinos from SN1987A have not decayed on their way 
one gets a lower bound on the electron neutrino lifetime as tq > 5.7 x 10^(mj,^/eV) 
sec. However if one includes neutrino mixing then shorter lifetimes are allowed 
provided |t/e2p < 0.81 [^. To be consistent with this, in our analysis we vary 
sin^ 9 in the range to 0.8. 
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The data used for the total rates is presented in Table 1. Since SK has much 
higher statistics, we have not included the Kamiokande results. The results ob- 
tained for Model 1 are summarised in Table 2. For the Model 2 z/2 decays to ui, 
which interacts as with a probabitity U^^ and with a probability U^2- We 
present in Table 3 the results for Model 2, where we have taken into account the 
i^e — e and — e scattering in addition to the t'e ~ e scattering in SK, while for the 
radiochemical experiments there is no change. Since the nui is degraded in energy 
m the effect is not significant. 



Experiment 


Chlorine 


Gallium 


Super-Kamiokande 


O bscrvcd Kate 
BP98 Prediction 


0.33 ±0.029 


0.562 ±0.043 


0.471 ±0.015 



Table 1: The ratio of the observed solar neutrino rates to the corre- 
sponding BP98 SSM predictions used in this analysis. The results are 
from refs. [0 and |l^. For Gallium, the weighted average of the SAGE 
and Gallex results has been used. 



Parameters 


Rates 
(d.o.f = 1) 


Spectrum 
(d.o.f = 15) 


Rates±Spectrum 
(d.o.f = 18) 


a (eV2) 





4.22 xlO-ii 


3.29 X 10-12 




0.5 


0.51 


0.32 


Xb 


Fixed at SSM value 


1.5 


0.72 


\^ ■ 

Amm 


12.59 


17.68 


33.59 


g.O.f 


3.88 xlO-2 % 


27.99% 


1.41% 



Table 2: The best-fit values of the parameters, the Xmin ^'^d ^^e good- 
ness of fit (g.O.f) for Model 1. 



Parameters 


Rates 
(d.o.f = 1) 


Spectrum 
(d.o.f = 15) 


Rates±Spectrum 
(d.o.f = 18) 


a (eV) 





3.3 xlQ-^' 


3.29 xlQ-''^ 


sin^e 


0.5 


0.53 


0.32 


Xb 


Fixed at SSM value 


1.5 


0.71 


A.mm 


11.98 


17.62 


33.12 


g.O.f 


5.38 xlO-2 % 


28.32% 


1.61% 



Table 3: The best-fit values of the parameters, the xLin the good- 
ness of fit (g.O.f) for Model 2. 
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From the analysis of only rates data the best-fit comes in the region where 
From the analysis of only rates data the best-fit comes in the region where a = 0, 
implying that the neutrinos are stable. For the total rates, the low energy neutrino 
flux should be suppressed less which is in contradiction with the energy dependence 
of the exponential decay term. Thus the best fit comes in the region where the 
decay term goes to 1. Since the probability in this region is energy independent 
the quality of the fit is not good. 

For the only spectrum analysis, the best-fit comes in the region where the 
exp(— term is non- vanishing and the high energy neutrinos are suppressed 
less. Since high energy bins have more number of events, the fit is much better 
compared to the fit to the total rates and is comparable to the ones obtained for 
the MSW oscillation solution (Xmm = 17.62) The vacuum oscillation solution 
gives a better fit 0. The best-fit values quoted in Table 2 are obtained with Xb 
constrained to be < 1.5. However if we remove the upper limit and allow Xb to 
take arbitrary values a slightly better fit is obtained for very high values of Xb'- 

• a = 6.53 X 10-" eV2, sin^ 9 = 0.8, Xb = 6.42, xLn = 17.22, g.o.f. = 30.41% 

We note that the a does not change much. But at higher Xb one needs a higher 6. 
To understand these features, in fig. 2 we plot Xb vs a for various fixed values of 
sin^ 9. For getting this curve we determine the Xb that gives the minimum at a 
particular a and then repeat this excersise for a varying in the range 10~^^ — 10^® 
eV^. The minimum obtained at each point of the parameter space of fig. 2 is 
within the 90% CL. limit of the global minimum. There are three regions 

1. For very low values of a (upto 10"^^ eV^) the exponential term is 1 and 
PueUe = 1 " 0-5 sin^26' which is the average oscillation probability. This can 
vary from 1 to 0.5 depending on 6. However in 11 out of 18 bins of the 
SK spectrum data, the rateobs/SSM < 0.5 and this can be achieved only by 
keeping Xb < 1-0. 

2. In the range 10^^'^ — 10^^° eV^ the exponential term contributes to the sur- 
vival probability and it falls sharply with increasing a. Therefore in this 
range Xb rises sharply as a increases. 

3. Beyond 10^^'^ eV^ the exponential term goes to zero and the probability is 
PueUe = cos'^6. In this zone one can achieve a probability < 0.5 by adjusting 
6 only and Xb does not play an important role. However if Xb is allowed to 
vary then for smaller values of 6 the Xb needed is low and vice-versa. 

In fig. 3 we plot the for the SK spectrum data, as a function of one of the 
parameters, keeping the other two unconstrained. In fig. 3a the solid(dashed) line 
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gives the variation of with a keeping Xb unconstrained (fixed at 1.0). If Xb 
is 1 then lower values of a are not allowed as by varying Q alone one cannot get 
a probability less than 0.5. If Xb is allowed to vary then low values of a are also 
admissible. Higher values of a are allowed for both the cases. The minimum of 
course comes in the region where a ~ 10~^^ eV^ for which the high energy neutrino 
events are less suppressed. 

Fig. 3b gives the variation of with sin^ Q keeping Xb unconstrained (solid 
line) and Xb = 1 (dashed line). For both the curves a can take any value. For 
Xb = 1 one cannot get a good fit in the low a region as discussed above and the 
fit goes to the high a region. In this zone the probability (~ cos^ 6) increases as 
we decrease 6 and this puts a lower limit on the allowed value of 6. This can be 
evaded if Xb is allowed to vary, as by adjusting Xb one can get a good fit even if 
6 is very low. 

Fig. 3c gives the allowed range of Xb keeping the other two parameters uncon- 
strained. From the figure we see that values of Xb below 0.4 are not allowed at 
90% C.L.. As discussed, Xb plays an important role in the low and intermediate 
a regions. In the latter zone the Xb required is high, therefore the constraints on 
the low values of Xb comes from the low a region. In this region as we decrease 
Xb the number of events will decrease which can be adjusted by increasing the 
survival probability and the maximum value of P^^ue = 1-0 gives a lower limit on 
Xb- For very high values of Xb the best-fit goes in the high a region where the 
probability is ~ cos'^6. As we increase Xb the number of events will increase which 
can be compensated by increasing 9. However from SN1987A constraints sin^ 9 is 
restricted to be < 0.8 and therefore beyond Xb = 7.2 one does not get a good fit. 

The rest frame lifetime obtained at the best-fit a from the spectrum data is 
To = (m2/eV)1.54 x 10~^ sec. This is small enough for decay to happen before 
neutrino decoupling in the early universe. This can increase the effective number 
of light neutrinos, N^, from 3. However depending on the data used the upper 



limit on N^, can be 5 or 6 ||I3[ which is consistent with the model of neutrino decay 
used here. 

The rest frame lifetime of 1^2 is given by [|l| 



167r 7712(1 + mi/m2) ^ 

^0 = — 7-^ (7) 

(7^ Am 



From the best-fit a ~ 10 eV^ and assuming 1712 » rrii one gets 

g^Am'^ ~ levrlO"" (8) 

If we now incorporate the bound (7^ < 4.5 x 10~^ as obtained from K decay modes 
T4| we obtain Am^ >~10~^ eV^, consistent with our assumption. At the best-fit 
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a the decay length for atmospheric neutrinos is A ~ 2.5 x 10^^ km and hence they 
do not decay. However for Am^ > 10~^ eV^ and large mixing angles there will be 
substantial v^. — conversion in conflict with the SK atmospheric neutrino and 



CHOOZ data 15, 16 



In conclusion, neutrino decay solution to the solar neutrino problem is ruled 
out at 99.96% CL. from the current data on total rates. For the SK spectrum data 
however, one can get much better fits, for lifetimes consistent with cosmological 
and supernova constraints. Although the best-fit for the spectrum data comes in 
the region where the probability is energy dependent, if Xb is allowed to vary 
the becomes flat over the entire range of parameter space. This implies that 
even energy independent suppression of the flux is allowed. Even if Xb is 
fixed at 1, the decay scenario can give comparable to the best-fit in the energy 
independent high a regime. For the flux normalisation factor a wide range 0.4 
< Xb < 7.2 is allowed at 90% CL. just from the SK spectral data, if neutrino decay 
is operative. This is much broader than the range allowed by SSM uncertainties. 
For the rate+spectrum analysis the decay solution is disfavoured at more than 
98% CL. even if Xb is allowed to take arbitrary values. 
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Figure Captions 



Fig. 1: The variation of the survival probabihty P^^i,^ (thick hues) and the ex- 
ponential decay term exp(— a L/E) (thin lines) with a, for two different values of 
neutrino energies. The solid curves correspond to neutrino energy = 7 MeV while 
the dotted curves are for neutrino energy = 13 MeV. For both the cases sin^ 6 is 
fixed at 0.51. 

Fig. 2: The variation of Xb with a for three different values of sin^^ shown in 
the plot. Each point on these curves is obtained by keeping a and sin^ fixed and 
determining the Xb corresponding to the minimum x^. 

Fig. 3: The variation of with (a) a for Xb unconstrained (solid line) and 
fixed at 1.0 (dashed line), while sin^6' is kept unconstrained for both curves; (b) 
sin^ 9 keeping both a and Xb unconstrained for the solid curve and with a un- 
constrained and Xb fixed at 1.0 for the dashed curve; (c) Xb keeping both a and 
sin^ 9 unconstrained. The dotted line shows the 90% CL. limit for 3 parameters 
(x^ = Xmtn + 6.25) and the dash-dotted line gives the corresponding limit for 2 
parameters (x^ = xLn + 4-61). 
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